Fuzzy recognition by the prokaryotic transcription factor HigA2 from Vibrio cholerae.

Disordered protein sequences can exhibit different binding modes, ranging from well-ordered folding-upon-binding to highly dynamic fuzzy binding. The primary function of the intrinsically disordered region of the antitoxin HigA2 from Vibrio cholerae is to neutralize HigB2 toxin through ultra-high-affinity folding-upon-binding interaction. Here, we show that the same intrinsically disordered region can also mediate fuzzy interactions with its operator DNA and, through interplay with the folded helix-turn-helix domain, regulates transcription from the higBA2 operon. NMR, SAXS, ITC and in vivo experiments converge towards a consistent picture where a specific set of residues in the intrinsically disordered region mediate electrostatic and hydrophobic interactions while "hovering" over the DNA operator. Sensitivity of the intrinsically disordered region to scrambling the sequence, position-specific contacts and absence of redundant, multivalent interactions, point towards a more specific type of fuzzy binding. Our work demonstrates how a bacterial regulator achieves dual functionality by utilizing two distinct interaction modes within the same disordered sequence.

Estimation of Peptide Helicity from Circular Dichroism Using the Ensemble Model.

An established method for the quantitation of the helix content in peptides using circular dichroism (CD) relies on the linear spectroscopic model. This model assumes an average value of the helix-length correction for all peptide conformers, irrespective of the length of the helical segment. Here we assess the validity of this approximation and introduce a more physically realistic ensemble-based analysis of the CD signal in which the length correction is assigned specifically to each ensemble conformer. We demonstrate that the linear model underestimates peptide helicity, with the difference depending on the ensemble composition. We developed a computer program that implements the ensemble model to estimate the peptide helicity. Using this model and the CD data set covering a broad range of helicities, we recalibrate CD baseline parameters and redetermine helix-coil parameters for the alanine-rich peptide. We show that the ensemble model leverages small differences in signal between conformers to extract more information from the experimental data, enabling the determination of several poorly defined quantities, such as the nucleation constant and heat capacity change associated with helix folding. Overall, the presented ensemble-based treatment of the CD signal, together with the recalibrated values of the spectroscopic baseline parameters, provides a coherent framework for the analysis of the peptide helix content.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins.

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Crystal Structure of de Novo Designed Coiled-Coil Protein Origami Triangle.

Coiled-coil protein origami (CCPO) uses modular coiled-coil building blocks and topological principles to design polyhedral structures distinct from those of natural globular proteins. While the CCPO strategy has proven successful in designing diverse protein topologies, no high-resolution structural information has been available about these novel protein folds. Here we report the crystal structure of a single-chain CCPO in the shape of a triangle. While neither cyclization nor the addition of nanobodies enabled crystallization, it was ultimately facilitated by the inclusion of a GCN2 homodimer. Triangle edges are formed by the orthogonal parallel coiled-coil dimers P1:P2, P3:P4, and GCN2 connected by short linkers. A triangle has a large central cavity and is additionally stabilized by side-chain interactions between neighboring segments at each vertex. The crystal lattice is densely packed and stabilized by a large number of contacts between triangles. Interestingly, the polypeptide chain folds into a trefoil-type protein knot topology, and AlphaFold2 fails to predict the correct fold. The structure validates the modular CC-based protein design strategy, providing molecular insight underlying CCPO stabilization and new opportunities for the design.

Structural basis of DNA binding by YdaT, a functional equivalent of the CII repressor in the cryptic prophage CP-933P from Escherichia coli O157:H7.

YdaT is a functional equivalent of the CII repressor in certain lambdoid phages and prophages. YdaT from the cryptic prophage CP-933P in the genome of Escherichia coli O157:H7 is functional as a DNA-binding protein and recognizes a 5'-TTGATTN6AATCAA-3' inverted repeat. The DNA-binding domain is a helix-turn-helix (HTH)-containing POU domain and is followed by a long α-helix (α6) that forms an antiparallel four-helix bundle, creating a tetramer. The loop between helix α2 and the recognition helix α3 in the HTH motif is unusually long compared with typical HTH motifs, and is highly variable in sequence and length within the YdaT family. The POU domains have a large degree of freedom to move relative to the helix bundle in the free structure, but their orientation becomes fixed upon DNA binding.

Binding-Induced Diversity of a Human Telomeric G-Quadruplex Stability Phase Space.

The structural polymorphism of G-quadruplex nucleic acids is an important factor in their recognition by proteins and small-molecule ligands. However, it is not clear why the binding of several ligands alters G-quadruplex topology. We addressed this question by following the (un)folding and binding of the human telomeric fragment 5'-(GGGTTA)3GGGT-3' (22GT) by calorimetry (DSC, ITC) and spectroscopy (CD). A thermodynamic analysis of the obtained data led to a detailed description of the topological phase space of stability (phase diagram) of 22GT and shows how it changes in the presence of a specific bisquinolinium ligand (360A). Various 1:1 and 2:1 ligand-quadruplex complexes were observed. With increasing temperature, the 1:1 complexes transformed into 2:1 complexes, which is attributed to the preferential binding of the ligand to the folding intermediates. Overall, the dissection of the thermodynamic parameters in combination with molecular modelling clarified the driving forces of the topological quadruplex transformations in a wide range of ligand concentrations and temperatures.

Analysis of Protein-DNA Interactions Using Isothermal Titration Calorimetry: Successes and Failures.

Isothermal titration calorimetry (ITC) is a golden standard for the characterization of protein-DNA binding affinities and allows direct assessment of the accompanying thermodynamic driving forces. Their interpretation can give insight into role of electrostatics, specificity of the DNA recognition, contribution of protein folding upon DNA binding and help to distinguish between minor and major groove binders. The main advantages of ITC are that the binding is measured in solution, and it requires no labeling of the samples, however, the method is not well suited for high-performance studies. Here we describe the sample preparation, a procedure to perform a typical ITC experiment, data analysis, and lastly discuss how to interpret the obtained thermodynamic parameters. In conclusion, we show examples of several unsuccessful ITC experiments and identify the underlying reasons for failed experiments. In most cases with a proper adjustment of the experimental setup, it was possible to obtain data appropriate for further analysis.

The free energy folding penalty accompanying binding of intrinsically disordered α-helical motifs.

Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes and preform critical roles in many cellular processes, most often through the association with globular proteins. Despite lacking a stable three-dimensional structure by themselves, they may acquire a defined conformation upon binding globular targets. The most common type of secondary structure acquired by these binding motifs entails formation of an α-helix. It has been hypothesized that such disorder-to-order transitions are associated with a significant free energy penalty due to IDP folding, which reduces the overall IDP-target affinity. However, the exact magnitude of IDP folding penalty in α-helical binding motifs has not been systematically estimated. Here, we report the folding penalty contributions for 30 IDPs undergoing folding-upon-binding and find that the average IDP folding penalty is +2.0 kcal/mol and ranges from 0.7 to 3.5 kcal/mol. We observe that the folding penalty scales approximately linearly with the change in IDP helicity upon binding, which provides a simple empirical way to estimate folding penalty. We analyze to what extent do pre-structuring and target-bound IDP dynamics (fuzziness) reduce the folding penalty and find that these effects combined, on average, reduce the folding cost by around half. Taken together, the presented analysis provides a quantitative basis for understanding the role of folding penalty in IDP-target interactions and introduces a method estimate this quantity. Estimation and reduction of IDP folding penalty may prove useful in the rational design of helix-stabilized inhibitors of IDP-target interactions. STATEMENT: The α-helical binding motifs are ubiquitous among the intrinsically disordered proteins (IDPs). Upon binding their targets, they undergo a disorder-to-order transition, which is accompanied by a significant folding penalty whose magnitude is generally not known. Here, we use recently developed statistical-thermodynamic model to estimate the folding penalties for 30 IDPs and clarify the roles of IDP pre-folding and bound-state dynamics in reducing the folding penalty.

Structural polymorphism of coiled-coils from the stalk domain of SARS-CoV-2 spike protein.

Spike trimer plays a key role in SARS-CoV-2 infection and vaccine development. It consists of a globular head and a flexible stalk domain that anchors the protein into the viral membrane. While the head domain has been extensively studied, the properties of the adjoining stalk are poorly understood. Here, we characterize the coiled-coil formation and thermodynamic stability of the stalk domain and its segments. We find that the N-terminal segment of the stalk does not form coiled-coils and remains disordered in solution. The C-terminal stalk segment forms a trimeric coiled-coil in solution, which becomes significantly stabilized in the context of the full-length stalk. Its crystal structure reveals a novel antiparallel tetramer coiled-coil with an unusual combination of a-d and e-a-d hydrophobic core packing. Structural analysis shows that a subset of hydrophobic residues stabilizes different coiled-coil structures: trimer, tetramer, and heterohexamer, underscoring a highly polymorphic nature of the SARS-CoV-2 stalk sequence.

Aggregation Time Machine: A Platform for the Prediction and Optimization of Long-Term Antibody Stability Using Short-Term Kinetic Analysis.

Monoclonal antibodies are the fastest growing class of therapeutics. However, aggregation limits their shelf life and can lead to adverse immune responses. Assessment and optimization of the long-term antibody stability are therefore key challenges in the biologic drug development. Here, we present a platform based on the analysis of temperature-dependent aggregation data that can dramatically shorten the assessment of the long-term aggregation stability and thus accelerate the optimization of antibody formulations. For a set of antibodies used in the therapeutic areas from oncology to rheumatology and osteoporosis, we obtain an accurate prediction of aggregate fractions for up to three years using the data obtained on a much shorter time scale. Significantly, the strategy combining kinetic and thermodynamic analysis not only contributes to a better understanding of the molecular mechanisms of antibody aggregation but has already proven to be very effective in the development and production of biological therapeutics.

Nanobody-aided crystallization of the transcription regulator PaaR2 from Escherichia coli O157:H7.

paaR2-paaA2-parE2 is a three-component toxin-antitoxin module found in prophage CP-993P of Escherichia coli O157:H7. Transcription regulation of this module occurs via the 123-amino-acid regulator PaaR2, which forms a large oligomeric structure. Despite appearing to be well folded, PaaR2 withstands crystallization, as does its N-terminal DNA-binding domain. Native mass spectrometry was used to screen for nanobodies that form a unique complex and stabilize the octameric structure of PaaR2. One such nanobody, Nb33, allowed crystallization of the protein. The resulting crystals belong to space group F432, with unit-cell parameter a = 317 Å, diffract to 4.0 Šresolution and are likely to contain four PaaR2 monomers and four nanobody monomers in the asymmetric unit. Crystals of two truncates containing the N-terminal helix-turn-helix domain also interact with Nb33, and the corresponding co-crystals diffracted to 1.6 and 1.75 Šresolution.

The sequence-ensemble relationship in fuzzy protein complexes.

Intrinsically disordered proteins (IDPs) interact with globular proteins through a variety of mechanisms, resulting in the structurally heterogeneous ensembles known as fuzzy complexes. While there exists a reasonable comprehension on how IDP sequence determines the unbound IDP ensemble, little is known about what shapes the structural characteristics of IDPs bound to their targets. Using a statistical thermodynamic model, we show that the target-bound ensembles are determined by a simple code that combines the IDP sequence and the distribution of IDP-target interaction hotspots. These two parameters define the conformational space of target-bound IDPs and rationalize the observed structural heterogeneity of fuzzy complexes. The presented model successfully reproduces the dynamical signatures of target-bound IDPs from the NMR relaxation experiments as well as the changes of interaction affinity and the IDP helicity induced by mutations. The model explains how the target-bound IDP ensemble adapts to mutations in order to achieve an optimal balance between conformational freedom and interaction energy. Taken together, the presented sequence-ensemble relationship of fuzzy complexes explains the different manifestations of IDP disorder in folding-upon-binding processes.

Unraveling the Thermodynamics of Ultra-tight Binding of Intrinsically Disordered Proteins.

Protein interactions mediated by the intrinsically disordered proteins (IDPs) are generally associated with lower affinities compared to those between globular proteins. Here, we characterize the association between the intrinsically disordered HigA2 antitoxin and its globular target HigB2 toxin from Vibrio cholerae using competition ITC experiments. We demonstrate that this interaction reaches one of the highest affinities reported for IDP-target systems (K D = 3 pM) and can be entirely attributed to a short, 20-residue-long interaction motif that folds into α-helix upon binding. We perform an experimentally based decomposition of the IDP-target association parameters into folding and binding contributions, which allows a direct comparison of the binding contribution with those from globular ultra-high affinity binders. We find that the HigA2-HigB2 interface is energy optimized to a similar extent as the interfaces of globular ultra-high affinity complexes, such as barnase-barstar. Evaluation of other ultra-high affinity IDP-target systems shows that a strategy based on entropy optimization can also achieve comparably high, picomolar affinities. Taken together, these examples show how IDP-target interactions achieve picomolar affinities either through enthalpy optimization (HigA2-HigB2), resembling the ultra-high affinity binding of globular proteins, or via bound-state fuzziness and entropy optimization (CcdA-CcdB, histone H1-prothymosin α).

Conditional Cooperativity in DNA Minor-Groove Recognition by Oligopeptides.

The recognition of specific DNA sequences in processes such as transcription is associated with a cooperative binding of proteins. Some transcription regulation mechanisms involve additional proteins that can influence the binding cooperativity by acting as corepressors or coactivators. In a conditional cooperativity mechanism, the same protein can induce binding cooperativity at one concentration and inhibit it at another. Here, we use calorimetric (ITC) and spectroscopic (UV, CD) experiments to show that such conditional cooperativity can also be achieved by the small DNA-directed oligopeptides distamycin and netropsin. Using a global thermodynamic analysis of the observed binding and (un)folding processes, we calculate the phase diagrams for this system, which show that distamycin binding cooperativity is more pronounced at lower temperatures and can be first induced and then reduced by increasing the netropsin or/and Na+ ion concentration. A molecular interpretation of this phenomenon is suggested.

A nanobody toolbox targeting dimeric coiled-coil modules for functionalization of designed protein origami structures.

Coiled-coil (CC) dimers are widely used in protein design because of their modularity and well-understood sequence-structure relationship. In CC protein origami design, a polypeptide chain is assembled from a defined sequence of CC building segments that determine the self-assembly of protein cages into polyhedral shapes, such as the tetrahedron, triangular prism, or four-sided pyramid. However, a targeted functionalization of the CC modules could significantly expand the versatility of protein origami scaffolds. Here, we describe a panel of single-chain camelid antibodies (nanobodies) directed against different CC modules of a de novo designed protein origami tetrahedron. We show that these nanobodies are able to recognize the same CC modules in different polyhedral contexts, such as isolated CC dimers, tetrahedra, triangular prisms, or trigonal bipyramids, thereby extending the ability to functionalize polyhedra with nanobodies in a desired stoichiometry. Crystal structures of five nanobody-CC complexes in combination with small-angle X-ray scattering show binding interactions between nanobodies and CC dimers forming the edges of a tetrahedron with the nanobody entering the tetrahedral cavity. Furthermore, we identified a pair of allosteric nanobodies in which the binding to the distant epitopes on the antiparallel homodimeric APH CC is coupled via a strong positive cooperativity. A toolbox of well-characterized nanobodies specific for CC modules provides a unique tool to target defined sites in the designed protein structures, thus opening numerous opportunities for the functionalization of CC protein origami polyhedra or CC-based bionanomaterials.

G-Quadruplex Stability from DSC Measurements.

Guanine-rich DNA oligonucleotides can adopt G-quadruplex (G4) structures in the presence of specific cations. Folding and unfolding of G4 can be characterized thermodynamically, providing the information on the stability of various G4 conformations. We show how the relevant thermodynamic and sometimes kinetic parameters are obtained by employing differential scanning calorimetry (DSC) and global fitting of an appropriate model to the DSC data. We demonstrate that best-fit values of the thermodynamic parameters can be interpreted in terms of the driving forces accompanying the G4 folding/interconversion and how they are translated into the phase diagrams, which provide an elegant description of the G4 phase space over a wide range of solution conditions.

A non-redundant data set of nanobody-antigen crystal structures.

A non-redundant data set of nanobody-antigen crystal structures is presented. The data set consists of a collection of cleaned pdb files which can be readily used as an input with most automatic analysis software. The accompanying data also include nanobody amino acid sequences with the annotated CDR regions. In the tabular format, we provide data on the interaction properties for each complex such as number of intermolecular interactions, experimental affinity and changes of the solvent accessible area. We also include the data regarding the surface composition of all nanobody and antigen molecules (surface occurrence of each amino acid type and its secondary structure). The data may be used for further structural bioinformatic studies of nanobodies and as the reference data when performing comparisons with the conventional antibodies.

Thermodynamic Stability of the Transcription Regulator PaaR2 from Escherichia coli O157:H7.

PaaR2 is a putative transcription regulator encoded by a three-component parDE-like toxin-antitoxin module from Escherichia coli O157:H7. Although this module's toxin, antitoxin, and toxin-antitoxin complex have been more thoroughly investigated, little remains known about its transcription regulator PaaR2. Using a wide range of biophysical techniques (circular dichroism spectroscopy, size-exclusion chromatography-multiangle laser light scattering, dynamic light scattering, small-angle x-ray scattering, and native mass spectrometry), we demonstrate that PaaR2 mainly consists of α-helices and displays a concentration-dependent octameric build-up in solution and that this octamer contains a global shape that is significantly nonspherical. Thermal unfolding of PaaR2 is reversible and displays several transitions, suggesting a complex unfolding mechanism. The unfolding data obtained from spectroscopic and calorimetric methods were combined into a unifying thermodynamic model, which suggests a five-state unfolding trajectory. Furthermore, the model allows the calculation of a stability phase diagram, which shows that, under physiological conditions, PaaR2 mainly exists as a dimer that can swiftly oligomerize into an octamer depending on local protein concentrations. These findings, based on a thorough biophysical and thermodynamic analysis of PaaR2, may provide important insights into biological function such as DNA binding and transcriptional regulation.

A dual role in regulation and toxicity for the disordered N-terminus of the toxin GraT.

Bacterial toxin-antitoxin (TA) modules are tightly regulated to maintain growth in favorable conditions or growth arrest during stress. A typical regulatory strategy involves the antitoxin binding and repressing its own promoter while the toxin often acts as a co-repressor. Here we show that Pseudomonas putida graTA-encoded antitoxin GraA and toxin GraT differ from other TA proteins in the sense that not the antitoxin but the toxin possesses a flexible region. GraA auto-represses the graTA promoter: two GraA dimers bind cooperatively at opposite sides of the operator sequence. Contrary to other TA modules, GraT is a de-repressor of the graTA promoter as its N-terminal disordered segment prevents the binding of the GraT2A2 complex to the operator. Removal of this region restores operator binding and abrogates Gr aT toxicity. GraTA represents a TA module where a flexible region in the toxin rather than in the antitoxin controls operon expression and toxin activity.

Energetic Basis of AGCGA-Rich DNA Folding into a Tetrahelical Structure.

It was recently discovered that, besides well-known G-quadruplexes and i-motifs, DNA may adopt another type of noncanonical structure called AGCGA-quadruplexes. Here, the folding of the VK2 fragment from the regulatory region of the PLEKHG3 gene is studied and, for the first time, the energetic contributions that stabilize this unique fold are described. Similarly to the B-DNA, it is stabilized by hydrophobic desolvation and, in contrast to G-quadruplexes, also by specific binding of water molecules. Compared to B-DNA, VK2 folding is enthalpically less favorable due to poorer base-stacking interactions, resulting in substantial conformational flexibility. This entropically favorable conformational "breathing" stabilizes the AGCGA-quadruplexes. In conclusion, AGCGA-quadruplexes have a distinguishing thermodynamic fingerprint and the corresponding driving forces enabling their folding are consistent with the observed structural features.